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Abstract 

We consider the Grand Unification (GU) scenario for neutrino masses which 
is based on the see-saw mechanism with the mass of the heaviest right handed 
(RH) neutrino at the GU-scale: M3 ~ Kqui and on the quark-lepton sym- 
metry for fermions from the third generation. The scenario predicts for the 
light neutrinos: 777,3 ~ (2 — 4) • 1CT 3 eV and m% ~ (0.3 — 3) • 10~ 5 eV (in the 
case of a linear mass hierarchy of the RH neutrinos or /and in presence of the 
Planck scale suppressed non-renormalizable operators). It also predicts large 
v e ~ v ii mixing: sin 2 29 e)1 > 0.2. In this scenario the solar neutrinos (vq) 
undergo both the v e — > v T resonance conversion in the Sun and substantial 
v e ~~ > v a vacuum oscillations on the way from the Sun to the Earth. The 
interplay of both effects enlarges the range of neutrino parameters which 
solve the z^-problem. In particular, v e — v T mixing angle can be as small 
as the corresponding quark mixing: sin 2 29 eT > (2 — 5) • 10~ 4 . The sce- 
nario predicts peculiar (oscillatory) distortion of the boron neutrino energy 
spectrum and seasonal variations of signals. Manifestations of these effects 
in the Super-Kamiokande and SNO experiments are studied. 



1. Introduction 

Small neutrino masses, in particular, the ones implied by the solar neutrino data, may be 
considered as a possible indication of Grand Unification (GU). This statement holds in the 
following context (which we will call the Grand Unification scenario for neutrino masses). 

1. Neutrino masses are generated by the see-saw mechanism 



MM' 1 j 

where w,d is the Dirac neutrino mass, = hoV u , hp is the neutrino Yukawa coupling, 

v u is the VEV of the Higgs doublet H u , and M is the Majorana mass of the right handed 

(RH) neutrino. A direct mass term for the left handed components may be generated by 

the induced VEV of the Higgs triplet of SU {2) L [0 . We assume that this term, if it exists, 

is much smaller than that in (JJ), at least for the heaviest neutrino. 

2. Quark-lepton symmetry is realized for the third generation of fermions at the Grand 

Unification scale, A GU = 2 • 10 16 GeV, which is defined as the scale of the gauge coupling 

unification in the Minimal Supersymmetric Standard Model. The quark-lepton symmetry 

relates the (Dirac type) Yukawa coupling of the neutrino with the Yukawa coupling of the 

top quark at Aqu'- 

h 3D (A G u) = h t (A GU ). (2) 

At the electroweak scale h t {MEw) — m t/v u , where m t is the top quark mass. The existence 
of quark-lepton symmetry and the equality (0) for the third generation are supported by the 
asymptotic mass relation m^Acu) — ^t{Aqu) which after renormalization effects results in 
a successful prediction of the 6-quark mass at low energies in supersymmetric GUT's ||. 

The boundary condition (|2|) is satisfied, in particular, in GUT's based on SOiq. Fur- 
thermore, in SOio-type unification all four Yukawa couplings in the third generation can be 
equal: h b = h T = h t = h 3 r>. 

Exact quark-lepton symmetry is, however, broken for the light generations. This follows 
from the failure of the asymptotic mass relations m s (Act/) = m^Acu) and m^A^t/) = 



m e (h.Gu) which are inconsistent with low energy determinations. One possible origin of this 
breakdown is non-renormalizable operators which are suppressed by Acu/Mp, where Mp is 
the Planck scale Q]. These operators can contribute significantly to the lighter generation 
masses while leaving the third generation relation essentially uncorrected. Therefore for the 
lighter generations one can expect relations of the type (|2[) to hold only as an order of mag- 
nitude. One would also expect the leptonic mixing angles involving the lighter generations 
to be different from their quark analogs. 

3. The Majorana mass of the RH neutrino from the third generation is at the GU 
scale: M 3 ~ Aqu- Using this mass, the boundary condition (§) and the known value of the 
top quark mass one finds from the see-saw formula the mass of the third neutrino (~ u T ) 
to 3 ~ several xl0~ 3 eV. The mass to 3 is in the range suggested by the Mikheyev-Smirnov- 
Wolfenstein (MSW) solution of the solar (z/ Q -) neutrino problem ||: 

to = (2 - 4) • 10~ 3 eV (3) 

(barring any degeneracy in the spectrum). Thus in the G£7-scenario the z/ Q -problem can be 
solved by the v e — > v T resonance conversion in the Sun provided the mixing between the 
first and the third generations is sufficiently large. Vice versa, using the value of the light 
neutrino mass (§) and the top quark mass at the electroweak scale along with the boundary 
condition @ one finds from the see-saw formula M 3 ~ 10 16 GeV ||. The approximate 
equality of the mass scales Aqu and M 3 may be considered as a hint for Grand Unification 
in addition to the gauge coupling unification and the b — r unification. 

An alternative solution to the solar neutrino problem is via v e — > resonance conversion. 
In this case the boundary condition for the second generation, m-iD^Gu) « to c (Ag{/) (to c 
is the charm quark mass), leads to M 2 ~ (10 10 — 10 11 ) GeV in the intermediate mass scale 
(see, e.g., 0). 

Note that in the 2z/-case, it is impossible to distinguish u e — > z/ M and u e — ► v T transitions 
using the z/ -data only. Both and v T are detected by the neutral current interactions 
which are the same (up to higher order corrections) for both neutrinos. The situation can 
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be different if the mixing of all three neutrinos is taken into account. We will show that 
solar neutrino data alone may disentangle the two possibilities (GU-scale scenario and the 
intermediate scale scenario). 

A possible discovery of //-oscillations by CHORUS/NOMAD || or new experiments like 
E803 (COSMOS) [[]] would imply m 3 ~ 0(1 eV), favoring the intermediate scale physics. 
Also confirmations of the LSND result and oscillation interpretation of the atmospheric 



neutrino anomaly will exclude the simplest version of the Gt/-scenario. 

In this paper we reconsider the Grand Unification scenario and study its signatures in 
the solar neutrinos. We calculate M 3 implied by the v e — > v T solution of the z/ -problem 
using the known mass of the top quark and discuss possible relations of M 3 to A GU (see sect. 
2). The v e - v T mixing is considered in sect. 3. We show that it is quite plausible that v e — > 
v T conversion will be accompanied by sizable v e — > z/ M vacuum oscillations of solar neutrinos 
(sect. 4). The interplay of both transitions (sect. 5) has a number of consequences: it 
enlarges the region of neutrino parameters in which one can get a correct description of the 
z/Q-data (sect. 6), it leads to peculiar distortion of the boron neutrino spectrum (sect. 7) 
and to seasonal variations of signals (sect. 8). Possible manifestations of these effects which 
can be considered as signatures of the Gt/-scenario in the Super-Kamiokande [ |TT| , |T2"| and 



SNO |13| experiments are studied in sects. 7, 8. In sect. 9 we summarize our main results 



and also comment on possible modifications of the GtZ-scenario which would allow one to 
explain other neutrino data. 

2. v e — v T conversion of solar neutrinos and the scale of Grand Unification 

Let us focus on the mass of the third neutrino m 3 which determines AmJ 3 responsible 
for the v e — v T conversion. We perform first a calculation of m 3 neglecting v T mixing with 
lighter generations. The procedure adopted is as follows. We use ([!]), (0) and the two-loop 
renormalization group equations (RGE) with the particle content of the Minimal Supersym- 
metric Standard Model. We fix the QCD gauge coupling a 3 (M z ) = g\{Mz)l^n = 0.118 at 
the Z° boson mass scale, M z , as an input and set the effective supersymmetry threshold at 
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M z . 



The renormalized neutrino mass at a scale \i < Kqu with the boundary condition (|2|) can 
be written as 



77l 3 (//) 



03 



(4) 



.k(A GI / 

where k is the coefficient of the effective dimension-5 neutrino mass operator, C e ff = 
KijLiLjH u H u , M03 is the mass of the RH neutrino in the absence of mixing between gener- 
ations. The RGE for k (which has been worked out only to one loop) is given by Jl4 



~dt 



167T 2 



6/i? + 2h 2 T - Q -g\ - Qg 2 2 



(5) 



Here g\ and g 2 are the U(l) and SU(2) gauge couplings, and h T is the Yukawa couplings of 
the tau lepton. The one-loop RGE for the top Yukawa coupling h t is 



dt 



ti 



6^ + ^ - Hp? - 3gl - fgl 



(6) 



where hf, is the bottom Yukawa coupling. 

In our numerical calculations of we kept the full two-loop RGE equations for h t , hb, 
h T . The results are summarized in Table I, where the values of m 3 are shown for different 
values of tan/3 = v u /v d and three possible values of m t ] the Majorana mass is fixed to be 
M 03 = 10 16 GeV. 

As follows from Table I the z/ r -mass is insensitive to the value of tan/5 except for tan/3 ^ 2. 
This feature can be understood from the one- loop semi-analytic expression for 7713. Indeed, 
the neutrino mass at a lower scale \x can be written as 



™?M ( a M \ A/m ( «3(X) 



m 3 (fx) 



-16/9 
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dt ( 3h 2 t + h\ 



hi 



(7) 



where o>g — 1/25 is the unified gauge coupling constant and a.\ = gf/4:7T. The neutrino 



mass depends on tan/5 via the Yukawa coupling h t = m t \Jl + tan 2 /?/(vtan/5). Obviously, 
for tan/5 > 2, the coupling h t is essentially independent of tan/?. 

Using the results of Table I and the mass of the light neutrino suggested by the solar 
neutrino data @ we can estimate the RH- neutrino mass M 03 . For m t = 175 GeV and 



tan/3 = 1.7 we find M 03 = (0.5 — 0.7)A G u- For tan/3 ^ 2 the top Yukawa coupling is near its 
fixed point value. Such a low value may be preferable from the point of view of successful 
b — t unification. On the other hand, for tan/3 > 2, we get M Q3 = (0.18 — Q.25)A G jj- Thus 
for the whole range of tan/? the Am 2 required for solution of the ^ Q -problem corresponds to 
M 3 in the interval 

M 03 < (0.2 - 0.7) A GU . (8) 

One remark is in order. If M03 is substantially below A G u, one should take into account 
the effect of the Dirac neutrino Yukawa coupling h 3 o on the evolution of h t , hb, h T and h 3 D 
itself in the interval between M 03 and A G u, since in this interval loops involving z/ 3 # will 



contribute [15|]. Although the momentum range for the running is rather small the effect 
can be substantial since both h t and h 3 r> are large (near the fixed point value of h t ). The 
boundary condition at M03 on the Dirac Yukawa coupling is modified to 



frjp(M 03 ) = hi D (A GU ) ( ai (M 03 ) \ 4/99 / 03(^03) r 16/9 

h 2 (M 03 ) ~ h 2 (A GU ) { a G ) \ a G ) 

1 r\nM 03 



X 



exp 



— dt (3h 2 3D - 3h 2 t +h 2 T -h 2 

Q1T Z JlnA GU K 



In the relatively small momentum range between the scales M03 and A G u the exponential 
factor above is almost unity because h 3 r> = h t and = h T at A G u. A difference between ht 
and h 3 o develops because of the difference in the gauge contribution. The effect of the gauge 
couplings is to diminish the Dirac mass m 3 E>, and consequently, the mass m 3 at low scale. 
Indeed, for M 03 /A G u = 1/10, m t = 175 GeV and the set of values tan/3 = (1.7, 3, 10, 30, 60) 
we find the neutrino mass m 3 = (2.15, 0.95, 0.82, 0.84, 0.96) • 10- 3 (10 16 GeV/M 3 ) eV. 
Comparing these values with the numbers in Table I we conclude that the threshold effect 
related to M 03 decreases the predicted value of m 3 by 5 - 25 %. Equivalently, the predicted 
value of M 3 for a fixed m 3 decreases by as much as 25% for tan/3 = 1.7 (fixed point value), 
and by 6% for moderate and large tan/? (M 03 /A G u = 1/10). So, the values of M Q3 are not 
changed significantly from those in (151). 
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Let us now take into account mixing of the third generation with the light generations 
of fermions. For simplicity we will consider mixing between two (say second and third) 
generations. Suppose Mo 2 is the Majorana mass of the RH neutrino component which 
leads via the see-saw mechanism to a certain mass m 2 in the absence of mixing between 
generations. For fixed values of m 2 and m 3 as well as the Dirac masses m 2 D, m 3 D, the 
mixing in the Majorana mass matrix of the RH neutrinos, M, will change values of the 
corresponding RH neutrino masses: M 2 7^ M02 and M3 7^ M03, where M 2 and M3 are the 
masses in presence of mixing. It can be shown, however, that the product of masses does 
not depend on the mixing, that is, M 2 ■ M3 = Mo 2 • M03 (up to small renormalization group 
corrections) J7j. From this equality we get a relation between hierarchies of masses in the 
presence of mixing, e = M 2 /M 3 , and without mixing, e = M 02 /M 03 : 

M 03 ^ 2 



The change M 03 — ► M 3 is related to the mixing angle 8m in M in the following way |7|: 

M 3 



sin 2 ^./ w e 



Mi 



(10) 



The mixing in M can raise the mass M3 to the GU - scale: M03 — > M3 ~ Agj/, for fixed 
masses of the light neutrinos and for fixed Dirac masses. This implies that the second mass, 
M 2 , decreases by factor M 03 /M 3 , and according to @ the hierarchy of masses is enhanced by 
a factor (M03/M3) 2 . If M 3 increases, e.g., by a factor 3, the hierarchy of masses is enhanced 
by one order of magnitude. Similar results can be obtained for mixing between the third 
and the first generations. 

From (8) we see that the mass of v 3 in the absence of mixing, M03, is (1.5 - 5) times 
smaller than Aqu- As follows from (|T0"| ) small mixing between generations 6m ~ 0(\/e) is 
enough to get the equality M 3 = Aqu. Another consequence of the increase in M 3 and the 
strengthening of the hierarchy is the enhancement of mixing of the light neutrinos which will 
be discussed in the sect. 3. 

The proximity of Aqu an d M 3 suggested by the solar neutrino data (or even possible 
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coincidence of these two values) can be considered as a hint for Grand Unification. Let us 
comment on possible relations between these two scales. 

Since the neutrino mass term Mv\vr is a singlet of the standard model symmetry group, 
there is no direct relation between the neutrino mass M and the scale of the gauge coupling 
unification. Some additional assumptions are needed to connect M 3 and Aqu an d, in general, 
one should not expect a coincidence of these scales. Possible relations between M 03 and 
Agu depend on the mechanism of Vr mass generation. Two simple mechanisms are worth 
mentioning: 

(i) The Majorana mass M3 is generated directly by the Yukawa couplings with a Higgs 
multiplet $ 

hsMVRVR® , (11) 

so that M3 = /i3m($)o- A relation between the interaction (|TT|) and physics of Grand 



Unification can arise in SOiq models, where <3> is the 126-plet. If $ is responsible for the 
symmetry breaking SO w — > SU 5 at the GU-scaAe, then ($) ~ A GU . Still h 3 M should be 
fixed. Note that h t (A GU ) ~ (0.4 — 1.5) for a wide range of tan/3. Therefore, to get M 03 in the 
region (^j) one should take h%M ~ (0.13 — 1.4) • ht in the absence of mixing. This allows for 
the possibility, especially if there is some mixing in the mass matrix M, that all the Yukawa 
couplings of fermions from the third generation are the same, in particular, h^M ~ h t . 
(ii) The mass M 3 can be generated by the effective operator 

^-/W^ > (12) 

where Ms is some mass at or above the G?7-scale. In this case M3 ~ h R (<&) 2 Q / Ms- In SO\q 
models $ is 16H-plet, and the operator (|T2| ) is generated by terms in the superpotential 

h R T6 H lQS + M S SS (13) 

which mix vr with an SOio-singlet or an adjoint, S, having the Majorana mass Ms- 
The 16^-plet breaks SO w to SU 5 and it is possible to identify its VEV with A GU . For 



Ms ~ (1 — 2) • Aqu one can take h^M ~ ht in agreement with equality of all Yukawa cou- 
plings in the third generation. 



3. v e — v T mixing 

In the two neutrino case the mixing angle needed to solve the z/ -problem should be in 
the range determined by [ TO : 

sin 2 29 = (3 - 10) • 1(T 3 . (14) 
As we will see in sect. 4, the effect of mixing between u e and v„ enlarges the region of 



solutions (0) to: 

sin 2 2# er = (0.2 - 10) ■ HT 3 . (15) 



Let us first compare (|T5|) with the corresponding mixing in the quark sector. Taking 9 er ~ 
V t d ~ (4 — 11) ■ 10~ 3 , where V t d is the element of the CKM quark mixing matrix, we get 
sin 2 29 eT = (0.06 — 0.5) • 10 -3 . These values cover the lower part of the range flr5|). Thus, 
the v e — v T mixing similar to the quark mixing is enough to solve the z/ - problem provided 
the v e — mixing is sufficiently large: sin 2 29 > 0.3 (see sect. 6). For sin 2 29 efl < 0.3 the 
values 9 er ~ Vtd are too small. In this connection let us consider the possibility of enhancing 
the mixing, so that 9 er ~ (2 — 3) Vu- 

In view of the violation of quark-lepton symmetry among lighter generations there is no 
reason to expect exact equality of 9 eT and V t d- A difference in quark and lepton mixing can 
come both from difference in the Dirac mass matrices: rri£> ^ m q (violation of the quark- 
lepton symmetry) and from the Majorana mass matrix (see-saw enhancement of lepton 
mixing |[7|])- Let us consider these two possibilities in order. 

1) The quark-lepton symmetry for light generations can be broken by the Yukawa cou- 
plings with additional Higgs multiplets. The contributions from these couplings to masses 
correct the bad asymptotic mass relations m^Acc/) = m e (AGu) and m s (Act/) = m^Acu)- 
Such corrections tend to enhance the mixing in the lepton sector. An explicit example has 
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been given in the context of SO\o ||18|| . There a minimal set of Higgs multiplets was intro- 



duced that couples to fermions: namely, a single 10-plet and a single 126-plet. The Standard 



Model singlet from the 126-plet generates the Majorana mass term for the RH neutrinos. 



The same 126-plet also contains a pair of SM doublets which receive vacuum expectation 
values, contributing to the quark and lepton masses. As a consequence, the structure of the 
Majorana mass matrix gets related to the quark and lepton mass matrices. The Yukawa 



couplings of 126-plet correct the bad asymptotic mass relations m^Acc/) = ^e{^Gu) and 
m s (AGu) — ^fiiA-Gu)- This led to the prediction 9 er ~ 3Vtd, where the factor 3 is the same 
Clebsch- Gordon coefficient that corrects the mass relations for light generations |19[]. In this 
case one gets sin 2 29 eT = (0.5 — 4) • - within the range fli~5|). 

2). An enhancement of lepton mixing can follow from mixing in the Majorana mass 
matrix of the RH neutrinos |I7|]. Indeed, the total lepton mixing angle (for two generation 
case) can be written as 

e = e D + e s , (16) 

where 9 b follows from the Dirac mass matrices of the neutrinos and charge leptons (9 = 9p, 
if M = I) and 9 S specifies the effect of the see-saw mechanism itself |T7]]. Let us consider 



mixing between the first and the third generation. The see-saw angle 9 S can be expressed in 
terms of the mass hierarchies = mi D /m 3D , e = Mi/M 3 and e = M i/M 03 [0: 



• 2 a ^ fc D 

sin u s ~ — 



/fo_ 1 



e 



M 3 1 



e LM 03 

For linear mass hierarchy of the RH neutrinos, eo ~ to, we get 



(17) 



sm 2 6 s « e D — 3 - . (18) 

Therefore €d = 10 -5 and M 3 /Mq 3 ~ 2 lead to 9 S ~ 3 • 10 -3 , which is comparable with Vta- 
If the hierarchy of masses is strengthened by one order of magnitude (e = O.leo), then the 
angle 9 S becomes 10~ 2 . In this case 9 S gives the main contribution to the lepton mixing, and 
sin 2 29 eT is of the order 10 -3 . 



According to fllPp for the linear hierarchy one gets 9m ~ \ft ~ \[^d- That is, the mixing 
in M similar to mixing in the Dirac mass matrices can lead to enhancement of the lepton 
mixing up to the value implied by the z/ -problem. Moreover, as we discussed in sect. 2, 
this enhancement can be related to increase of M3 to A.qu- 

There is an alternative mechanism for generating enhanced lepton mixing. Mixings in 
the lepton sector can differ substantially from those in the quark sector as a consequence of 
the nature of Grand Unification symmetry breaking. The stronger hierarchy observed in the 
up-quark masses relative to the hierarchies in the down quarks and charged lepton masses 
gives some indication towards such a possibility. One example was suggested in |20 |. The 



main idea is that in SU 5 (or in the SU 5 decomposition of SO w ) Ui and ul quarks lie in the 
same 10j representation, while the di and d\ (and similarly leptons l\ and U ) are in separate 
representations: 10j and 5j respectively. If a hierarchy factor is associated with the 10- 
plets, but not with the 5-plets, the down quark and charged lepton masses (5 -10 ) will 
scale as e^, whereas the up-quark masses (10 -10) will scale as e 2 D . The left-handed CKM 
angles will be small, as they scale as eo. The right-handed CKM mixing angles for quarks 
(unphysical in the Standard Model) will be of order one, since the 5i have no hierarchy 
factor. Since the left-handed lepton doublets are in 5,, their mixing angles will be large. 
The specific model constructed in Ref. leads to 9 eT ~ 0.032 which is within the range 



4. Parameters of v e — system 

If the mixing angles in the Majorana matrix M are small, the mass of the second neutrino 
can be estimated as 

where h c is the Yukawa coupling of the charm quark. For M02 ~ M03 ~ Aqu we find 
m,2 < 10~ 7 eV and Amj 2 w 10 -14 eV 2 which is far below the sensitivity of the z/ -data. 
Let us assume a linear mass hierarchy in the RH neutrino mass matrix, according to which 
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the eigenvalues of M have the same hierarchy as the eigenvalues of the Dirac mass matrix 
m D : Mi oc m^, or e ~ e . In this case M 02 ~ 10 14 GeV, and for the light neutrino we get 
from the see-saw formula (P m 2 ~ (0.3 — 3) ■ 10~ 5 eV. Consequently, 

Am 2 12 « (10- 11 - 1(T 9 ) eV 2 , (20) 

which is in the range where vacuum oscillations on the way from the Sun to the Earth are 
important. 

Let us comment on the possible origin of the linear mass hierarchy of M. If in SOio 



models the same 126-plet contributes both to the RH neutrino masses and to the charged 



fermion masses, the hierarchy is naturally linear [£L[ . 

The linear hierarchy may be due to family symmetry. Let us consider a U(l) symmetry 
whose breaking is characterized by a single small parameter A (which may be the VEV of 
a singlet Higgs scalar divided by the Planck mass) having the [/(l)-charge —1. We assume 
that the fermionic 16j-plets of 50(10) carry family U(l) numbers gj. Then as a consequence 
of the U(l) invariance, the mass term 16jl6j will be suppressed by a factor \ qi+q 3 . Since 
the Majorana mass terms ((|11|) or (|T2"D) have the same flavor structure as the Dirac mass 
terms they will have similar hierarchy in the eigenvalues provided that the only source of 
£/(l)-breaking is A. 

Note that in the case of linear mass hierarchy the lightest RH neutrino has the mass 
(10~ 6 — 10~ 5 ) Acu ~ (10 10 — 10 11 ) GeV which is in the correct range to explain baryogenesis 
via leptogenesis p2| . 

The u e — v ^ mixing angle is expected to be of the order the Cabibbo angle, 9 C . It can 
be larger than 8 C due to violation of exact quark-lepton symmetry in the light generations. 
Furthermore, an additional enhancement may follow from the see-saw mechanism itself, as 
was discussed in sect. 3. Therefore 9 efM ~ (1 — 2) 9 C , and consequently, sin 2 2# eM ~ 0.2 — 0.7 
are quite plausible without any additional assumptions. This mixing can lead to observable 
effects in the solar neutrinos. 

In the G£7-scenario the masses of v\ and v 2 are so small even in the case of the linear mass 
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hierarchy of the RH components that possible contributions from the Planck scale physics 
can be important. In general, an expression for neutrino masses is the sum 

m « m s + m nr , (21) 

where m s is the see-saw mass and m nr is the contribution from possible effective non- 
renormalizable interactions associated to the Planck scale physics [ [23)1 : 



LjLiH u H u . (22) 



M P 

Here ctij are constants expected to be of order 1, are the leptonic doublets and Mp is the 



Planck mass. The interaction ( P2|) may follow immediately from string compactification or 
from renormalizable interactions via the exchange of particles with mass ~ Mp. 

The interaction (|22|) gives the mass m nr ~ ~TJ~ v u ~ 10~ 5 eV. The contribution m nr 

IVlp 

to the mass of the third neutrino can be neglected, whereas m nr is important for and 
it dominates in m\. Again, Am 2 2 ~ m 2 r ~ 1CT 10 eV 2 is in the range fl2Tj) . Also the 
mixing is modified. The interactions ( f2"2] ) determine the mixing between the first and second 
neutrinos. If a, L j ~ 0(1), this mixing angle can be large: sin 2 29 e ^ ~ 0(1). The mixing 
of light generations with the third generation remains small: e.g., 9 eT ~ m nr /m^ ~ 0.01. 
Nevertheless this contribution is comparable to V t d- Therefore the total v e — v T mixing can 
be enhanced by interactions (0), so that sin 2 26 er will be in the interval fllip. 



5. The interplay of resonance conversion and vacuum oscillations 

As we have established in sects. 3, 4 the GU- scenario with linear mass hierarchy of 
the RH neutrino masses or/and with additional effects from the non-renormalizable Planck 
scale interactions can naturally provide the pattern of the neutrino masses and mixing with 
m 3 ~ (2 - 3) ■ 10~ 3 eV, m 2 ~ (0.3 - 3) • 10" 5 eV, m 1 < m 2 , sin 2 29 eT = (0.2 - 3) ■ 10~ 3 
and sin 2 29 efl ~ 0.2 — 0.7. For these values of parameters both the v e — > v T resonance flavor 
conversion and v e — vacuum oscillations on the way from the Sun to the Earth become 
important. Moreover, non-trivial interplay of these two effects takes place. We will call such 
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a possibility the hybrid (resonance conversion plus vacuum oscillation) solution of the solar 
neutrino problem. 

Due to the mass hierarchy and the smallness of 9 er the dynamics of the three neutrino 
system is reduced to a two neutrino task. Indeed, the electron neutrino state can be written 
in terms of the mass eigenstates Vi (i = 1, 2, 3) as 

v e = cos9 er -i> + sin# eT -z/ 3 , (23) 

where 

v = cos 6 efl ui + sin # eM z/ 2 . 

Inside the Sun due to the smallness of Am^, the system v\ — v-i is "frozen". The evolution 
of the whole system consists of the v e resonance conversion to the state v' = cos9 eT us — 
sin 9 eT v (orthogonal to v e ). On the way from the surface of the Sun to the Earth the state 
z/3 decouples from the system: large mass difference Am^ leads to averaged oscillation effect 
or/and to loss of coherence. On the way to the Earth the v e — vacuum oscillations occur 
due to mass splitting Am 2 2 between v\ and u 2 . Taking this into account it is easy to write 
the v e survival probability pi : 



P = P v (Ami 2 , g-P^Am^, 6 eT ) + 0(sm 2 6 eT ). (24) 

Here Py(Am^ 2 , 9 efi ) is the 2v - vacuum oscillation probability characterized by parameters 
Am^ 2 , 9 efl , and Pr is the 2v - averaged survival probability of the resonance conversion 
characterized by Amf 3 and 9 er (Pr is averaged over the production region). For sin 2 2# eT < 
3- 10~ 3 the 0(sin 2 9 eT ) - corrections in (^) can be safely neglected and the total probability is 
factorized. The exact formula has been derived in . In fig. 1 we show typical dependence 
of the survival probability on the neutrino energy. 

The survival probability ( p4]) satisfies the following inequalities: 

(1 - sin 2 29 etl ) -P R < P < P R . (25) 

That is, P is an oscillatory function of the neutrino energy inscribed in the band fl2"5| ) (fig. 
1). The width of the band equals AP = sin 2 29 efJi ■ Pr. 
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If vacuum oscillations are averaged out, we get from (|2 

P = (1 - ^sm 2 26 efl )-P R . (26) 
The properties ( [25] , p6| ) allow one to derive immediately several consequences of the 



interplay of the vacuum oscillations and resonance conversion 



The solar v e are transformed both to and v T : the corresponding transition proba- 
bilities equal P{y e — > v T ) « 1 — P R and P{y e — > f M ) ~ Pr — P, where P and Pr are 
defined in (1241). 



• According to (p5|) , vacuum oscillations lead to additional suppression of the v e -flux 
in comparison with the pure resonance conversion. As a consequence, new regions 
of neutrino parameters appear in which one can describe the i/ -data. In particular, 
sin 2 29 eT < 1CT 3 are allowed now (sect. 6, fig. 2). 

• For Am 2 2 < 1(T 9 eV 2 the vacuum oscillations are not averaged and lead to additional 
oscillatory distortion of the 8 B- or/and pp-neutrino spectra (sect. 7). 

• One expects additional seasonal variations of neutrino flux due to the dependence of 
the vacuum oscillation probability on the distance from the Sun to the Earth (sect. 
8). 

• For large Am 2 2 the pp- neutrino flux is suppressed by the averaged vacuum oscillation 
effect. For small Am 2 2 < 1CT 11 eV 2 the suppression depends on the neutrino energy, 
thus leading to distortion of the pp-neutrino energy spectrum. 

Thus, a detailed study of the energy spectra of the 8 B- neutrinos, and (in future) of pp- 
neutrinos as well as measurements of seasonal variations of the signals will allow one to test 
the Gf/-scenario. 

6. New regions of parameters 
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The interplay of the resonance conversion and vacuum oscillations opens new possibilities 



in description of the i/ -data |24j]. Some of these possibilities were discussed in [25|| . In 
particular, for sin 2 29 eT ~ 10~ 3 and Am^ ~ 10 -4 eV 2 the high energy part of the boron 
neutrino spectrum lies in the nonadiabatic edge of the two neutrino suppression pit (due 
to the resonance conversion). For these values of parameters the resonance conversion does 
not change neutrino fluxes at low and intermediate energies. The suppression of the 7 Be- 
neutrino flux implied by the z/ -data can be due to vacuum oscillations, if e.g., the first high 
energy dip of the oscillation probability is at E ~ 1 MeV. For this to occur, Am 2 2 should 
be in the range (7 - 9)-10~ 12 eV 2 . Thus, the vacuum oscillation suppression pit is at low 
energies, and the resonance conversion pit is at high energies. Note that Am 2 3 ~ 10 -4 eV 2 
correspond to values of M3 substantially below the GU-scaXe: M3 ~ 5 • 10 14 GeV. 

In contrast, for the Gt/-scenario with Am 2 3 = (4 — 10) • 1CT 6 eV 2 the 7 Be -flux can 
be suppressed by the resonance conversion. Since sin 2 29 eT ^ 1CT 3 , the suppression pit is 
narrow and the high energy (observable) part of the boron neutrino spectrum is suppressed 
insufficiently An additional suppression can be due to vacuum oscillations, especially if the 
observable part of the boron neutrino flux is situated in the first high energy dip of the 
oscillatory curve (see fig. 1). Thus we get the configuration with the resonance conversion 
pit at low energies and the vacuum oscillation pit at high energies. The region of parameters, 
where this configuration gives a good fit of the data is shown in fig. 2. (In the x 2 -analysis we 
considered the Super-Kamiokande and Kamiokande as separate experiments. Thus, there is 
one degree of freedom for the hybrid solution and three degrees of freedom for 2z/ solution 
of the solar neutrino problem.) As follows from fig. 2 values of mixing angles as small as 
sin 2 29 eT = 2 • 10~ 4 are now allowed. 

Further diminishing of 6 er implies an increase of 8 efJi in such a way that the hybrid 
solution converges to pure vacuum oscillation solution. Furthermore, it is impossible to 
diminish Am 2 3 and therefore to increase M3 in the hybrid solution. The minimal value, 
Am 2 3 « 4 • 10~ 6 eV 2 , is defined by the condition that 7 5e-neutrino line is at the adiabatic 
edge of the resonance conversion pit. 
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For the pp-neutrinos one gets the averaged oscillation effect: P ~ 1 — 0.5sin 2 2# eM ~ 
0.7 — 0.9. Therefore, the Germanium production rate in the Gallium experiments is typically 
at the lower border of the allowed region (~ 60 SNU's). This suppression is weaker than the 
one in the pure vacuum oscillation solution. 

Since the mixing angle responsible for the resonance conversion is very small, sin 2 29 er < 
3 • 10~ 3 , the day-night effect is negligible. 

7. Distortion of the boron neutrino spectrum and signals in 
SuperKamiokande and SNO 

An interplay of vacuum oscillations and resonance conversion can lead to peculiar dis- 
tortion of the boron neutrino energy spectrum. In particular, for Am 2 2 > 10 -11 eV 2 one 
expects an additional oscillatory modulation of the spectrum due to vacuum oscillations (fig. 



Let us consider a manifestation of such an oscillatory distortion in the energy spectrum 
of the recoil electrons measured by the Super-Kamiokande as well as (in future) SNO exper- 
iments. 

In presence of the neutrino conversion described by the survival probability P{E V ) the 
number of recoil electrons, N(E vis ), with a given visible energy, E vis , equals: 



where E e is the total energy of the recoil electron, <&(E V ) is the original boron neutrino flux, 
f(E vis ,E e ) is the energy resolution function which can be parameterized as 



1). 




P{Eu) 



da Ve 
dE e 



(E e , E v ) + (1 - P(E v ))-^(E e , E v ) 



(27) 




(28) 



We use the value of a determined in the Super-Kamiokande calibration experiment [[n|,|T2|. 
It has an approximate dependence on the neutrino energy a oc -7==. 
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Similar number of events No(E V i S ) = N(E V i S ,P = 1) has been found in the absence of 
conversion. Using N and N we calculate ratios, R l e , of events with and without oscillations 
in the AE = 0.5 MeV energy bins corresponding to the Super-Kamiokande presentation: 

™ = fi! +AE dE 'vis N ( E 'vis) (9Q , 

e IS +AE dE> is N (E' vis ) • { 3) 

In figs. 3, 4 we show (by histograms) the ratios R l e expected for different values of the 

neutrino parameters as well as the ratios measured by the Super-Kamiokande experiment 

during 306 days ||12|| . For comparison we present also the distortions expected from the pure 

MSW-solution (fig. 3a) and pure vacuum oscillation solution (fig. 3b). 

As follows from fig. 4, the integrations over the neutrino energy and the electron energy 
convoluted with the resolution function lead to strong averaging of the oscillatory behavior. 
Indeed, for present water Cherenkov detectors, like the Super-Kamiokande, the width of the 
resolution function (on semi-height), 2a, is ~ 4 MeV at 10 MeV which is comparable with 
the width of the largest dip of the oscillatory curve (in the energy scale). To illustrate the 
effect of smoothing we show the ratios for the ideal resolution f(E V i S ,E e ) = 5(E v i s — E e ) . 
The most profound effect follows from the first (the widest) dip of the oscillatory curve. 

Let us consider manifestations of the oscillatory behavior in details. For this we fix 
parameters of the v e — v T system and find the modification of the dependence of R e on E vis 
for different values AmJ 2 . We will fix sin 2 29 efl = 0.5. An increase of the angle 9 e)1 leads to 
enhancement of distortion. 

With increase of Am\ 2 the oscillatory curve which corresponds to the vacuum oscillation 
probability shifts to high energies (see fig. 1). In such a way different parts of the oscillatory 
curve turn out to be in the observable range of the boron neutrino (z^b-) spectrum: E ~ 
(5 - 15) MeV. 

We find the following distortion picture for different intervals of Am 2 : 

• Am 2 2 < 10~ n eV 2 : The oscillatory curve is below the observable range (5 — 15) MeV. 
Therefore the distortion of ^-spectrum is close to that produced by the resonance 
conversion alone. 
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Am 2 2 < 6 • 10~ n eV 2 : the i/B-spectrum is on the rising (with energy) part of the first 
(high energy) dip of the oscillatory curve. In this case the vacuum oscillations enhance 
the distortion produced by the resonance conversion. The slope of the R e (E V i S ) is 
bigger than in the pure MSW-case (fig. 4 a). 

AmJ 2 = (6 — 10) ■ 1CT 11 eV 2 : the minimum of the first dip of the oscillatory curve is 
in the range (5 — 15) MeV. This leads to the dip in the R e (E vis ) dependence (fig. 4 b, 
c). 

Am 2 2 = (1.0 — 1.2) • 10~ 10 eV 2 : the z/^-spectrum is in the decreasing part of the first 
dip which manifests itself as the negative slope in the distortion of R e (E V i S ) (fig. 4d). 

Am 2 2 = (1.2 — 2) • 10~ 10 eV 2 : the first (high energy) bump of the oscillatory curve is 
in the range (5 — 15) MeV which produces the bump in the recoil electron spectrum 
(fig. 4e, 4f). 

Am 2 2 = (2 — 5) ■ 10~ 10 eV 2 : the second dip of the oscillatory curve is in the range 
(5 — 15) MeV. Again with increase of Am 2 2 first an enhancement of the slope occurs 
(fig. 4g), then a dip appears, and finally the slope becomes negative (fig. 4h). 

Am 2 2 > 5 ■ 10 -10 eV 2 : the second dip shifts to E > 15 MeV. There are several narrow 
dips in the observable part of spectrum. Some non-averaged effect still exists in the 
high energy part. 

Am 2 2 > 10~ 9 eV 2 : The periods of the oscillatory curve in the observable interval are 
substantially smaller than 2a. Strong averaging of the oscillation effect takes place. 
The ratio R e approaches the asymptotic dependence which can be obtained from (|27D 
by substitution Py — > Py = 1 — 0.5 sin 2 29 efM : 

R e ^R?-P v + (l- Pv)^r = (1 - 0.5 sin 2 29) ■ R? + 0.5 sin 2 29^ , (30) 

where Rf is the ratio for the case of pure resonance conversion, is the effect of 
muon (tau)-neutrino: 
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= J dE e J dEJ-®-^. (31) 

Note that this contribution does not depend on the probabilities of transitions. Nq 
is the number of events in absence of oscillations. The last term in (|3H| ) is relatively 
small: ~ 0.04. 

According to ( |3~0| ) for large Am^, vacuum oscillations result in flattening of the dis- 
tortion stipulated by the resonance conversion alone (decrease of the slope). 

In certain intervals of Am^ an additional effect of vacuum oscillations leads to distortion 
which contradicts observations. Therefore, in these intervals large mixing angles (typically 
sin 2 29 efM > 0.5) can be already excluded. 

The smoothing effect is weaker in the SNO experiment. The integration over the neutrino 
energy gives weaker averaging and energy resolution is expected to be slightly better. In 
our calculation we use the cross-sections of reaction vd — > epp from |2B[ and the energy 



resolution function in the form (28) with a(E e ) = 14%/y 10MeV/-E e - One needs experi- 
ments with at least two times better energy resolution (like HELLAZ |27|) to measure the 
oscillatory modulation of the energy spectrum. 

Deviations of the spectrum distortion from the simple form predicted by the pure MSW 
solution or vacuum oscillation solution will indicate the effect of the third neutrino considered 
above. 

8. Seasonal variations of signals 

The vacuum oscillation probability depends explicitly on the distance between the Sun 
and the Earth L: Py = Py(L). Therefore seasonal variations due to pure geometrical effect, 
1/L 2 , related to the eccentricity of the Earth orbit will be modified [28"|. Depending on 



values of the neutrino parameters, the geometrical effect can be enhanced or suppressed by 
the L-dependence in the probability. Moreover, for experiments which are sensitive to the 
high energy part of the boron neutrino spectrum the modification is strongly correlated to 



the distortion of the recoil electron energy spectrum |29 
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Let us consider the seasonal variations in the case of hybrid solution. Using the proba- 
bility ( f24|) we can rewrite the expression for the number of events in the Super-Kamiokande 
experiment (|27| ) in the following way: 

da Ue da Vll 



N 



J dE e J dE u f-$-P R -P v - 



+ iV M , (32) 



dE e dE e 

where is defined in ([[!]). Note that iV M changes with time only due to the geometrical 
dependence of the flux. If the resonance conversion probability varies weakly in the observ- 
able energy interval, we can substitute it by certain average value Pr. Then the expression 
for the number of event becomes 

NnP R - N vac + (1 - P R ) ■ , (33) 

where N vac is the number of events expected in the case of pure vacuum oscillations. Notice 
that for the configuration shown in the fig.l the probability Pr indeed changes in a small 
interval from 0.7 to 0.9. According to (|33| ) the vacuum oscillation effect is attenuated by the 
factor Pr. This means that to get the same overall suppression of flux one needs smaller 
value of sin 2 29 efl , and therefore the seasonal variations become weaker. Moreover, there is 



an additional damping term in 

Let us introduce the seasonal asymmetry A as 

N w -N s 

= 2 NW TW > (34) 

where N w = J w dtN(t) and N s = J s dtN(t) are the integral numbers of events during the 
winter and summer time correspondingly. Using (|33|) we get 

N„ 



A = A 



v ■ 



i - a - ft) ■ # 



(35) 



Here Ay is the asymmetry in the case of pure vacuum oscillations and iV is the number of 
events averaged over the year. The last term in ( pop is typically smaller than 10%. As we 
mentioned above, in presence of the resonance conversion the value sin 2 29 eil which leads to 
a good fit of the data is smaller than in the pure vacuum oscillation case: sin 2 2# eM < 0.8. 
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Then from (|35|) we find that seasonal variations related to the probability are attenuated by 
some factor < 0.7 in comparison with variations in the pure vacuum oscillations solution. 

9. Discussion and conclusions 

Let us first summarize our main results. 

1. We have considered the Grand Unification scenario which is based on the see-saw 
mechanism with the Majorana mass of the RH neutrino M3 ~ Aqu and on the quark-lepton 
symmetry for fermions from the third generation. The Gt/-scenario leads to the mass of 
the third neutrino, m 3 ~ (2 — 4) • 10 -3 eV, in the range of masses implied by the resonance 
conversion of the solar neutrinos. 

2. We have found that in the G£7-scenario the MSW - solution of the ^ -problem implies 
the mass M 3 ~ (0.2 — 0.7)Aqu (depending on the value of tan/3) in the absence of mixing 
of the RH neutrinos and M 3 can coincide with Aqu if small {6 m ~ ^ M1/M3) mixing in the 
Majorana mass matrix of the RH neutrinos is taken into account. 

The proximity of the mass M3 to Aqu can be considered as an indication of the Grand 
Unification. 

3. Assuming a linear hierarchy of the RH neutrino masses or /and the presence of 
contributions from non-renormalizable operators related to Planck scale physics we get 
m 2 ~ (0.3 — 3) 10~ 5 eV 2 and Am 2 2 ~ (10~ n — 10~ 9 ) eV 2 . Furthermore, the mixing between 
the first and the second generations can be rather large: sin 2 26 efJL > 0.2. In this case the 
v e — > vacuum oscillations on the way from the Sun to the Earth give rise to observable 
effects. 

4. The G^-scenario leads to non-trivial interplay of the v e — > v T resonance conversion 
and v e —>■ vacuum oscillations of the solar neutrinos. An additional vacuum oscillation 
effect enlarges possible range of neutrino mixing, so that sin 2 26 er can be as small as 2 • 10~ 4 
which is of the order of corresponding quark mixing. 

5. The interplay of the resonance conversion and vacuum oscillations leads to new type 
of distortion of the boron neutrino spectrum, namely, to oscillatory modulation of the energy 
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dependence produced by the resonance conversion alone. 

The integrations over the neutrino energy and the electron energy convoluted with the 
energy resolution function in the Super-Kamiokande experiment (and similarly SNO) result 
in strong smoothing of the oscillatory behavior of the recoil electron energy spectrum. Still 
depending on Am^ 2 one can observe an enhancement of the positive slope, or appearance of 
negative slope as well as dip or bump in the dependence of R e on E vis . Two times better 
energy resolution is needed to observe real oscillatory behaviour. 

Certain ranges of parameters are already excluded by existing data. 

6. Seasonal variations of signals are expected due to dependence of the vacuum oscillation 
probability on distance from the Sun to the Earth. Typically these variations are expected 
to be weaker than in the pure vacuum oscillation case. 

The GU- scenario considered in this paper does not allow one to explain the LSND result 
or to solve the atmospheric neutrino problem in terms of neutrino oscillations. It predicts 
"null" results in the CHORUS and NOMAD experiments as well as in future short and 
long baseline oscillation experiments. Null result is also expected in the neutrinoless double 
beta decay searches etc. Therefore a positive result in at least one of these experiments 
will exclude the simplest version of the G£/-scenario. Also the scenario does not supply any 
appreciable contribution to the hot dark matter (HDM) in the Universe. 

It is possible to modify the G£/-scenario discussed here so that the second neutrino 
mass is m 2 ~ (0.1 — 1) eV. For this the mass of the corresponding RH neutrino should be 
M 02 ~ 10 10 GeV in the absence of mixing (or for small mixing). Such a modification can 
accommodate the LSND result and give a sizable contribution to HDM. An explanation of 
the atmospheric neutrino anomaly requires also large mixing between and u T . This is 
more difficult to achieve, since strong mass hierarchy accompanied by large mixing angle 
implies further tuning of parameters. 

Additional freedom in explaining of the data arises if new (sterile) neutrino states are 
introduced. 
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TABLES 

TABLE I. The mass of third neutrino 7713 as the function of tan (3 for M R = 10 16 GeV and 



three possible 


values of physical top quark mass: 


m t = 175, 170, 180 GeV. 




tan/3 




v T -mass (10~ 3 eV) 






m t = 170GeV 


m t = 175GeV 


m t = 180GeV 


1.7 


1.47 


2.85 


— 


1.9 


1.16 


1.64 


2.. 71 


3 


0.83 


0.99 


1.17 


5 


0.76 


0.88 


1.00 


10 


0.73 


0.84 


0.95 


30 


0.74 


0.85 


0.96 


50 


0.76 


0.92 


1.05 


60 


0.83 


0.99 


1.17 
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FIGURES 




FIG. 1. The averaged (over production area) survival probability P{y e — ► u e ) (solid line) 
as a function of the neutrino energy E for the following values of parameters: sin 2 29 eil = 0.77, 
Am 2 .! = 6.0 • 10- n eV 2 , sin 2 26 er = 8.0 • 10~ 4 , Ara^ = 1.1 • 10~ 5 eV 2 . The dependence of the 
averaged MSW probability PR{v e — > v e ) alone on E is also shown (dashed line). 
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95% CI. 




FIG. 2. The 95% C.L. regions of solution of the z/ -problem in Am^ and sin 2 26 eT plot (306 
days Super-Kamiokande data are taken into account) for different vacuum oscillation parameters 
sin 2 26 L, Am 2 !, and the original boron neutrino flux (in units of the standard solar model) 
= <Z>( 8 B)/$ssm( S B). The thin solid lines correspond to sin 2 26 efl = 0.48, Am^ = 1.4 • 10" 10 eV 2 , 
and the values of /b indicated at the curves. The thick lines correspond to (a) sin 2 29 efl = 0.66, 
Am^ = 8.3 • 10~ n eV 2 and / B =1.0; (b) the same as (a) but /b=0.9; (c) sin 2 20 e ^ = 0.62, 
Am 2 .! = 9.5 • 10- 11 eV 2 , f B = 0.8; (d) sin 2 29 efl = 0.50, Am|j = 7.8 • 10^ n eV 2 , f B = 0.7. 
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FIG. 3. The expected spectrum deformations of the recoil electrons in the Super-Kamiokande 
(solid lines) and SNO (long dashed lines) experiments for (3a) two neutrino conversion with 
sin 2 20 = 8.8 • 1(T 3 , Am 2 = 5.0 • 10" 6 eV 2 , (3b) vacuum oscillation with sin 2 29 = 0.82, 
Am 2 = 6.4 • lO" 11 eV 2 . The curves are normalized so that i? e (10MeV)=0.368 for the Su- 
per-Kamiokande and i? e (10MeV)=0.256 for the SNO experiment. Thin lines correspond to ideal 
energy resolution; thick lines correspond to real energy resolution. The Super-Kamiokande 306 
days data are shown (statistical errors only). 
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Amii =3.0x1 0" 11 eV 2 



■ SK: detector E resolution 
SK: ideal E resolution 
SNO: detector E resolution 
SNO: ideal E resolution 



9 10 11 12 13 14 15 
Evis (MeV) 

(4a) 



Anrii =6.0x1 0" 11 eV 2 



0.70 - 
0.65 
0.60 
0.55 
0.50 
0.45 - 
0.40 - 
0.35 
0.30 - 
0.25 - 

0.20 rr. 

0.15 — 



SK: detector E resolution 

SK: ideal E resolution 

SNO: detector E resolution 

SNO: ideal E resolution 




10 



(4b) 
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Amii =8.0x1 0" 11 eV 2 
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Amii =6.0x10 10 eV 2 
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FIG. 4. The same as in figs. 3 but for the hybrid solution with parameters sin 2 29 efM = 0.5, 



sin 2 29 eT = 6.0 • 10 4 , Am^ = 8.0 • 10 6 eV 2 and different values of Am 2 i indicated in the figures. 
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